Exact Kohn-Sham eigenstates versus quasi-particles in simple models of strongly correlated 
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We present analytic expressions for the exact density functional and Kohn-Sham Hamiltonian of simple tight- 
binding models of correlated electrons. These are the single- and double-site versions of the Anderson, Hubbard 
and spinless fermion models. The exact exchange and correlation potentials keep the full non-local dependence 
on electron occupations. The analytic expressions allow to compare the Kohn-Sham eigenstates of exact density 
functional theory with the many-body quasi-particle states of these correlated-electron systems. The exact Kohn- 
Sham spectrum describes correctly many of the non-trivial features of the many-body quasi-particle spectrum, 
as for example the precursors of the Kondo peak. However, we find that some pieces of the quasi-particle 
spectrum are missing because the many-body phase-space for electron and hole excitations is richer 
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I. INTRODUCTION 

Density Functional Theorjiii^ allows to tackle complex 
quantum systems comprising N interacting electrons. Its 
essence consists on the replacement of the extremely convo- 
luted many-particle electronic interactions with an effective 
one-body potential, also known as the exchange and correla- 
tion potential V^'~'' , whereby the A^-particle Hamiltonian is 
substituted by a fictitious one-particle Kohn-Sham Hamilto- 
nian H^^ . The exact V-^'-^ is however not known and it is 
a widespread belief that it is not possible to find an analytic 
exact expression for it. The popularity of Density Functional 
Theory (DFT) has arisen from the fact that semi-empirical fit- 
tings of V^'~^ to the exchange and correlation potential of the 
Jellium model in the Local Density Approximation (LDA) and 
improvements over it^^- perform remarkably well for a large 
majority of materials, molecules and nanostructures. The 
qualitative features of the structural and functional properties 
of many systems are usually well reproduced, and in a number 
of cases, quantitative agreement can also be reached. How- 
ever, these practical implementations of DFT are not perfect. 
They fail to predict a number of relevant properties, specially 
for strongly correlated electronic systems. Yang and cowork- 
ers have discussed some explicit conditions that exact energy 
functionals must obey^. 

DFT has also been proposed for tight-binding models of 
strongly correlated electrons^.. The availability of exact semi- 
analytical or numerical results for the ground state energy 
as a function of the electron concentration in the Hubbard 
and the spinless fermion models^i^ has allowed to establish 
a Bethe ansatz LDA theory for the m^'"i'^ . An extension of 
the theory to describe time-dependent external potentials has 
enabled the description of non-equilibrium electron transport 
phenomenai^— . However, Bethe ansatz LDA theory also has 
limitations. First, since the Bethe Ansatz solution expresses 
the ground-state energy in terms of the electron concentra- 
tion, only a local density approximation could be formulated. 
Second, the analytic formula for ground state energy is exact 



only at half-filling, while away from it semi-analytical or nu- 
merical fittings to the solution of the Bethe Ansatz equations 
must be performed. Finally, inhomogeneous systems where 
strong correlations take place in localized region of them are 
better described by Anderson models. However, the Ander- 
son model is Bethe Ansatz-solvable only if the band of un- 
correlated electrons is linearizedi^, leaving the system energy 
unbounded from below. Therefore, the ground state energy 
can not be obtained by a minimization procedure, which ren- 
ders the Bethe ansatz LDA approach useless for the Anderson 
model. 

The quasi-particle (QP) excitation spectrum of a system de- 
termines its response to external perturbations according to 
Landau's Fermi liquid theory. Furthermore, this spectrum 
is directly accessible via spectroscopic techniques of differ- 
ent sorts. It would therefore be quite useful if the Kohn- 
Sham (KS) eigenstates provided at least a qualitative descrip- 
tion of it, as one would expect to happen at least for systems 
where electronic correlations are weak. This is indeed con- 
firmed by a vast amount of calculations and comparisons be- 
tween KS eigenvalues and experimental or numerical data of 
weakly coiTelated materials. However, quantitative agreement 
is sometimes not so good. Furthermore the KS spectrum is 
frequently qualitatively wrong in strongly-correlated materi- 
als. Notice now that even if the exact V-^'-^ of a specific sys- 
tem is known, a possible correspondence between the exact 
KS and the exact many-body QP spectra is not supported at 
all by the basic theorems of DFT. An exception is the Highest 
Occupied Molecular Orbital (HOMO) which by Janak's the- 
orem is given by the chemical potential of the system, which 
is a ground-state property^—. In other words, DFT predicts 
the correct position of the HOMO level of a system, provided 
that the exact V-^'-'' , or a very good approximation to it, is 
known. Failures to predict the correct position of the HOMO 
must therefore be attributed to a poor approximation to the ex- 
act V-^'-' . Failures to reproduce the rest of the spectrum could 
however be due either to limitations of DFT proper, or to a 
poorly approximated functional. Indeed, while an exact func- 
tional may not provide a good description of the full QP spec- 
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trum, it is clear that if in addition, the quality of the approxi- 
mate functional is poor, the proposed spectrum of KS eigen- 
values will bear a small resemblance to the true QP spectrum. 
Since no exact functional for a strongly correlated system has 
ever been developed, the above two sources of disagreement 
have never been fully disentangled. The main goal of this arti- 
cle is to separate them. We will find the exact KS eigenvalues 
of several simple models of strongly correlated electrons and 
compare them with the exact many-body QP spectrum. This 
will allow us to understand the size of the self-energy correc- 
tions to the exchange-correlation potential. 

One of the main sources of disagreement between approx- 
imate DFT KS eigenvalues and exact QP originates in the 
mean-field-like treatment of electronic correlations which lie 
at the heart of LDA. Indeed, electrons behave as quantum 
point particles. However, mean-field theories replace quan- 
tum probabilities by classical clouds of charge. As a result, 
every electron may interact with its own charge cloud, giv- 
ing rise to spurious direct and exchange self-interaction ef- 
fects. Additionally, each electron interacts with the clouds of 
other electrons having opposite spin, leading to what is some- 
times called the static correlation error^. For systems contain- 
ing more than one atom, these mean field clouds are spread 
throughout the whole entity in contrast to electrons which 
are always point particles and therefore fully localized. This 
spread gives rise to further spurious effects termed delocaliza- 
tion errors, which lead to incorrect dissociation energies and 
QP excitation energies for molecules^. A prototypical exam- 
ple of the delocalization error is an iJ^ molecule in the dis- 
socation limit where the two ions are held widely apart. The 
single electron in the molecule has equal probability of resid- 
ing in any of the two atoms, but a measuring proccess will 
find it fully localized in only one of them. Mean field theories 
in contrast place half an electron in each ion. The excitation 
energy of an added quasi-electron will therefore be different 
in the two cases— '^'. 

Improving the description of the QP spectrum therefore im- 
plies improving the description of electronic correlations. The 
Hartree-Fock approximation as well as the self-interactions 
correction schemo^i22r21 get rid of the self-interaction effect, 
but not of other mean-field drawbacks. To go beyond these 
schemes, the Dyson-Sham-Schliiter equation must be used 
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where G^p^^^ is the Green's function obtained from the ap- 
proximate KS Hamiltonian. Notice that G^pp,.^.^ carries al- 
ready a mean-field description of the electron interaction. A 
perturbative expansion for the self-energy E"^'-^ must then 
be set up to improve the description of correlations and 
in particular to amend the destruction of quantum effects 
brought about by the mean-field approximation. The GW 
approximation^*" has been quite successful in the description 
of electronic and optical properties directly linked to the QP 
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but does not correct the problems mentioned 



above. Some recent work by Romaniello and coworkers show 
how the careful inclusion of vertex corrections allows to get 
rid not only of the self- interaction effects but also of part of the 
delocalization effects22i212^. However, delocalization debris 



remains since molecular dissociation is still not well handled. 
In addition, Millis and coworkers22i2£ have studied the per- 
formance of the GW approximation for the Anderson model, 
and shown how this approximation can describe Coulomb 
blockade effects, but fails to describe the emergence of Kondo 
Physic a^'i^^ . Dynamical Mean Field Theory, implemented to- 
gether with an accurate impurity-solver, includes many of the 
most relevant short-range correlation effects^"—. 

We have devised a procedure that has allowed us to find an- 
alytic expressions for the exact energy density functional of 
the single- and double-site Anderson, Hubbard and spinless 
fermion models, from which we have been able to write down 
the corresponding exact Hamiltonians H^^. Since the QP 
spectrum of these models is available analytically from con- 
ventional many-body techniques, we have been able to per- 
form explicit and detailed comparisons of the full spectra of 
exact KS eigenvalues and of exact many-body QP. We have 
found that the KS eigenvalue corresponding to the HOMO 
level agrees with the corresponding QP state. This implies that 
the exact H^^ correctly predicts that the lowest energy for 
electron addition of an A^-electron system is equal to the high- 
est energy for electron removal of the corresponding A^ + 1 
system^*'. We have also found that the exact H^^ of the An- 
derson model describes correctly the emergence of the Kondo 
resonance and of other quasi-particles. However, we find that 
there exact density functional theory misses pieces of the ex- 
act many-body QP spectrum. A way to improve the descrip- 
tion of the spectrum would be to use again the Dyson-Sham- 
Schliiter equation 
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where G^^^^j is the Green's function associated to the ex- 
act KS Hamiltonian. We expect that this self-energy and its 
perturbative expansion should be much simpler than the self- 
energy defined in Eq.(l) above, because now the unperturbed 
Green's function retains the full quantum nature of electrons. 
Our piece of work is complementary to efforts by other groups 
to provide exact functionals for simplified systems. We men- 
tion here recent work by Burke and collaborators, who have 
found numerically exact density functionals for some one- 
dimensional models by combining DFT with Density Matrix 
Renormalization Group techniques^. 

The layout of this article is as follows. Section II describes 
the methodology employed to find out exact functionals for 
systems with a small number of electrons. This methodol- 
ogy is applied in sections III and IV to describe the single-site 
Hubbard model, and the double-site Anderson model, respec- 
tively. The conclusions are laid down in section V. The solu- 
tion of the double-site Hubbard model is placed in appendix 
A. The solution of the double-site spinless fermion model can 
be found in appendix C. 



II. METHODOLOGY 

We begin with a description of our method, which is based 
on the formulation of DFT on a lattice^. We have found 



that the conventional ensemble-based method to describe non- 
integer occupationsl^i^ fails in the formulation of the exact 
density functional of the single-site model described below. 
We have therefore devised an alternative method which is 
specifically adapted for the description of quantum systems 
with a small but not necessarily integer number of electrons 
N. 

We consider a physical system whose time-evolution is dic- 
tated by a tight-binding Hamiltonian. As an example, we write 
explicitly the Hamiltonian of the Anderson model. 






where a set of N electrons hop back and forth along a chain 
of i = 1, .., A^ atoms, labeled by the index c, and to another 
atom, denoted by the index d, where electron correlations take 
place via a Coulomb term U. The cr— index denotes the up and 
down components of the electron spin. 

We use the Fock space of states of the system {\(j) > } to set 
up our variational scheme. Site occupations, electron numbers 
and the expectation value of the Hamiltonian are given by 
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We do not use a Hartree term in the definition of E^^ be- 
cause we have found no traces of such a term in the analytic 
equations for the exact functionals. Therefore we have found 
it useless for the purposes of the present discussion. We define 
, the exact KS Hamiltonian as follows: 
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We wish to define an energy density functional 
Q[nc,i,cr,nd,a',U] whosc minimization gives the exact 
ground state energy E'^ and occupations n°g. for a target 
set of electron numbers {NS,N?). To define Q, we note 
that every given set of occupations {nc,i,cT,nd,a'} can be 
reproduced by several states \(f> >. In other words, if we 
classify these states in boxes labeled by each occupation set, 
then each box contains several | <f> >, and each of these has 
a different energy E {(/)). However, if we choose in each box 
the state | (/)™ > with minimum energy E"^ = E{(j>™-), we 
achieve a one-to-one correspondence between occupation 
sets and energies for every box, which allows to define the 
energy density functional (5[nc,i,cr, J^d,;!', C^] = i?™~. Since 
there exist in general several sets of occupation numbers 
{?^c,i,cr,?^d.cr'} giving the same target electron numbers 
N^ = iV°, the ground state energy E'^ is obtained by 
minimizing Q over all those sets. This procedure then defines 



the ground state occupations {n° ^ 
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We define now the non-interacting kinetic energy functional 
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i,U = 0], and the Exchange- 



correlation functional E-^'-^ = Q — T, from which the ex- 
act Exchange-correlation potential V^'~' is obtained by taking 
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where Edc is a double-counting term. Notice that the above 
procedure allows to define functionals and KS Hamiltonians 
for systems with a fractional electron number However, the 
many-body Hamiltonian in Eq. (3) commutes with the elec- 
tron number operator A^^- Therefore the many -body Hamilto- 
nian eigenstates must describe an integer number of electrons, 
unless some degeneracy occurs. We will see later on that the 
functional Q has a polygonal shape, so that the Exchange- 
correlation potentials jump by constants at integer N^ values, 
which lead to ambiguous definitions of the KS eigenvalues 
at integer A^^r- However, the total energies of the ground- 
and excited-states of the KS Hamiltonian E" {a = 0, 1, ...) 
are continuous because the the jumps in the summations over 
KS eigenvalues are counterbalanced by similar jumps in the 
double-counting terms. The ground state energy _E" of the ex- 
act KS Hamiltonian and many-body Hamiltonians agree with 
each other by construction, but this is not so for the total en- 
ergies of the excited states of both Hamiltonians, which are 
needed to construct the Green's functions. 

The QP spectrum of the many-body Hamiltonian can be 
compared with the KS and mean-field spectra by looking 
at the poles and residues of the Green's functions G{u!), 
G^^^{ll!) and G^'^^{uj). We define on this matter the many- 
body HOMO level as the QP peak which is partially filled. 
The exact G and G^^ need not agree, except for the pole 
describing the HOMO level. We will use the Lehmann 
representation"*'^ to compute G for integer A^- values. In ad- 
dition, the equations-of-motion method"" yields a closed set 
of equations for G for the single-site model. This method 
nicely enables to extrapolate the G-poles to non-integer elec- 
tron numbers, and agrees with the results obtained using the 
Lehmann representation for integer A^. The mean field spec- 
trum can be obtained from the eigenvalues of the one-body 
mean-field Hamiltonian, or using the equations-of-motion 
method for G^^^ . The KS spectrum could also be obtained 
from the eigenvalues of the one-body KS Hamiltonian. How- 
ever, these KS eigenvalues are discontinuous at A^ integer so 
ambiguities in the ascription of eigenvalues to QPs arise for 
integer A^ . It is therefore essential to use the Lehmann repre- 
sentation as a guide. 



We close this section by describing an alternative proce- 
dure which also allows us to find exact results. If the ex- 
act ground state energy E'^ and occupations n'^^ are found 
then the Schrodinger equation for the Kohn-Sham hamilto- 
nian can be inverted to find the exact exchange and correla- 
tion potential corresponding to the ground state occupations 



,,xc 



K- 



XClmO 



We note however that vf^^ is not a func- 
tional, but rather corresponds to the Exchange-correlation po- 
tential functional evaluated at the ground state occupations. 
This procedure is simpler than the methodology described in 
this section, but does not allow us to find functionals. Similar 
methods have been employed by Baerend and coworkers'*-, as 
well as by Helbig and coworkersi^ to find exact analytical or 
numerical expressions for the Exchange-correlation potential 
of diatomic molecules in the dissociation limit. 



III. SINGLE-SITE ANDERSON-HUBBARD MODEL 
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FIG. 1: (Color online) Three-dimensional plot of the exact energy 
functional Q of the single-site Anderson-Hubbard model as a func- 
tion of (wd.t) '^d.i), for a value of ej = 1 and of C/ = 10 (in arbitrary 
units). 



The above methodology can be easily applied to the single- 
site Anderson-Hubbard model (TM = 0), where only two oc- 
cupations {rid.tj^^d,^} are defined. The states in the Fock 
space of the single-site model can be expressed using the num- 
ber basis |nd,-f,"-d,; > as 

I > == ao I 0, > + a^ I 1, > + a^ | 0, 1 > + 02 | 1, 1 > 

(7) 
The expectation value of the occupation numbers and the en- 
ergy can then be expressed as 
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where D ~ |aop + |a^p + |a^p + |a2p. The simplest way 
to find Q[nd.^,nd,i] is as follows. Solve first for some of 
the coefficients a^ using the occupations for the occupation 
numbers n^ cr- Those coefficients are then eliminated by in- 
serting them back into the equation for E{(t>). The resulting 
expression is then minimized in terms of the remaining coef- 
ficients. One must be careful though to choose coefficients 
which are strictly non-zero in a given domain of N . In the 
present case, it is best to solve for a^ since these are finite for 
all A^ = Ud.^ + 7id.i different from 0, 2: 

(1 - N) (|atP + la^n = N |ao|' + (2 - N) |a2p (9) 



The resulting equation for 
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is minimized as 
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The above expressions for the exact functional Q can be 
summarized as 



Qind.t^ "dj = edN + U{N-1)0{N- 1) 



(12) 



where 9 is the Heaviside step function. This expression gives 
the correct ground state energy for a system with a target num- 



ber iVJ', A^[ of electrons 
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Note that this ground state energy is spin-degenerate. Sub- 
tracting from Q the non-interacting kinetic energy functional 
T[nd,-f,nd,i] = edN, and taking a functional derivative, 
we find the exact Exchange-correlation potential V^J^ = 
U 9{N - 1). The resulting KS Hamiltonian 



iJ^^^ = 



Z^ 
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provides the correct E^ thanks to the double counting term 
U6{N° - 1). Notice that the KS eigenvalue jumps by U 
exactly at A^ = 1, and is therefore ill-defined at that inte- 
ger A^-valu Q-^i'^-'^'^"^ . The density functional Q has the cor- 
rect trapezoidal shaped as a function of No-, from which the 
right expression for the chemical potential of the system can 
be obtained. Furthermore, Q shows flat-plane behavior when 
plotted as a function of the occupation numbers, as displayed 
in Fig. 1. We note that Yang and coworkers established some 
exact conditions on the shape of the exact energy functional, 
from which they deduced such a flat-plane behavior^. These 
conditions enabled them to draw an educated plot of the en- 
ergy functional of the hydrogen atom, which is very similar to 
our Fig. 1. 

We write down now the mean-field Hamiltonian of this 
model 
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FIG. 2: (Color online) Quasiparticle spectrum of the single-site 
Anderson-Hubbard model for M = 0. The shaded gray area shows 
the position of the many-body poles, where the area is proportional 
to the weight of the peak. The black solid line represents the location 
of the exact Kohn-Sham eigenstate. The poles of the paramagnetic 
mean-field solution are shown with a dashed red line. Energy units 
are arbitrary. 



where we have subtracted the conventional mean-field double- 
counting term. The mean-field Hamiltonian gives the follow- 
ing estimate for the energy of the system 
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where the spin-degeneracy of the exact solution is lost. Note 
that in the Hubbard and Anderson models every electron inter- 
acts only with electrons of opposite spin. As a consequence, 
the mean-field theory does not suffer from direct or exchange 
self-interaction effects. However, because of the mean-field 
replacement of electron probabilities by charge clouds, an 
electron of spin a interacts with a fraction Ud.^a of electrons 
of opposite spin, instead of with a full electron of opposite 
spin with probability Ud.-a- As a consequence, in the para- 
magnetic solution M = 0, every electron interacts artificially 
with a fraction N/2 of electrons of opposite spin. However, 
the mean-field ground state energy is minimized by the fully 
spin polarized solutions M = N, in which case the spurious 
interaction between opposite-spin charge clouds is avoided by 
a wrong mechanism and, as a consequence, E^^^ ^ E^ . In 
contrast, the interacting piece of the exact KS Hamiltonian 
U 6{N — 1) is only activated if a full electron exists already in 
the system, and therefore retains the full quantum behavior. 

The exact many-body and mean-field QP spectrum are ob- 
tained from the poles and weights of the retarded Green's 
function 
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The above equations can easily be obtained using the 
equations-of-motion method and allow to extrapolate the QP 



spectrum to non-integer values of nd,cr which, coupled to the 
spin degeneracy of the total energy enable the exploration of 
different spin states. 

We analyze first the paramagnetic state where M = and 
nd.a = N/2. The exact KS Green's function is found by 
combining the equations-of-motion method with the Lehmann 
representation for TV = 0, 1, 2. The following formula sum- 
marizes the results 
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and extrapolates them to non-integer- iV-values. The many- 
body, exact DFT and mean field Green's function are shown 
in Fig. 2 as a function of N . The many-body Green's func- 
tion has two poles, which can be viewed as the ancestors of the 
lower and upper Hubbard bands of the Hubbard and Anderson 
models^. The position of these two poles depends neither on 
the occupation nor on the spin of the system. They are sepa- 
rated exactly by an energy U and their weight shifts smoothly 
from one peak to the other as N increases. Because exact 
DFT is a single-particle theory, its Green's function yields a 
single peak per KS eigenvalues, whose weight equals one. A 
remarkable exception happens at iV = 1, where the KS eigen- 
value show an abrupt change from e^ to ed + U ■ Notice that 
both eigenvalues contribute at A^ = 1 with equal weight. The 
A^ ^^ 1^ (A^ -^ 1+) KS eigenvalue exactly agrees with the 
many-body lower (upper) Hubbard band precursor. This non- 
trivial result allows to draw an important conclusion: even if 
KS eigenvalues show a jump at integer electron number val- 
ues, the eigenvalues at both sides of the given A^ contribute to 
the QP spectrum. To summarize, the positions and weights of 
the exact KS and many-body peaks coincide for integer num- 
bers A^, showing how exact DFT keeps the quantum nature of 
the electrons in spite of being a one-body theory. The abrupt 
shift at A^ = 1 can be viewed as the way that exact DFT uses 
to retain that quantum nature: if there is less than one electron 
at the site, then there is no Coulomb interaction because the 
electrons does not interact with itself. If there is more than one 
electron, then the Coulomb interaction between point particles 
of opposite spin is activated, rising the energy by U . Notice 
that many -body and exact DFT agree on the value of HOMO 
level, which is also equal to the chemical potential [i, defined 
as the derivative of the total energy with respect to the particle 
numberi^. 

We must remember however that in this quantum system 
only states with integer electron numbers A^^r = 0, 1 are mean- 
ingful. Therefore for A^ = 1, the system must contain a full 
electron with either spin up or down. We therefore turn now 
to analyze a maximally spin-up polarized. The up- and down- 
spin Green's functions are different now, as is apparent from 
the Lehmann representation at A^ = 1, where we take 1 1, > 
as the ground state. 
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FIG. 4: (Color online) The eight pieces in the (A'^, A''_i,)-plane which 
must be used to perform the constrained minimization procedure 
leading to the exact Q-functional of the double-site Anderson and 
Hubbard models. The thick dots indicate the positions where G and 
G^ are evaluated. 



FIG. 3: (Color online) Quasi-particle spectrum of the single-site 
Anderson-Hubbard model for a maximally spin-up polarized case, 
(a) and (b) show the poles of the spin-up and spin-down Green's 
functions, respectively. The shaded gray area represents the many- 
body quasi-particle, where the width is proportional to the peak 
weight. The black thick dot represents the location of the exact Kohn- 
Sham eigenstates. The poles of the spin-polarized mean-field solu- 
tion are shown with a dashed red line. Energy units are arbitrary. 



where Ejq denote the total ground state energy for N = 
0, 1, 2. We determine now the many -body and mean field QP 
spectra for fractional occupation numbers using Eqs. (16). To 
determine G^^"^ correctly for integer N we use the Lehmann 
representation. The following formula extrapolates G^^ to 
non-integer A^-values 



G^f(^) 



G^f(^) 



9{1- N + S) 
UJ ~ ed + iS 

9{N -l + S) 
UJ- {ed + U)+iS 



(20) 



The different spectra are shown in Fig. 3. As before, the exact 
KS eigenvalues agree with the many-body QP for integer A^. 
The mean-field states have a closer resemblance to the many- 
body QP, although clear differences still exist, whose origin 
is traced back to the static correlation error As a closing re- 
mark, we note that the GW approximation cures these mean- 
field artifacts for the present case as shown by Romaniello and 
coworkersaii. 



IV. DOUBLE-SITE ANDERSON MODEL 

The model in the previous section has allowed to show 
how the exact density functional retains the quantum nature 



of electrons in a single atom and therefore avoids the static 
correlation error brought about by mean field theory. We wish 
to address in this section how the exact functional avoids also 
the delocalization error in a strongly-correlated model con- 
taining two sites. We show that the exact KS Hamiltonian 
provides a correct description of the atomic limit of the cor- 
related model. The discussion is centered in the Anderson 
model, but our conclusions can be also applied to the double- 
site Hubbard model, which is solved in Appendix A. Notice 
that the many-body QP spectrum of the full Anderson model 
is much more complex than that of the single-site model dis- 
cussed above and, in addition to the lower and upper Hubbard 
bands it develops a Kondo resonance in the Kondo regime. 
We therefore wish to explore now whether exact DFT could 
describe this more convoluted QP spectra. Finally, notice that 
the obtained energy functional Q is spin-degenerate, in con- 
trast to the full Anderson and Hubbard models, where this de- 
generacy is absent. It is therefore interesting to check whether 
exact DFT lifts the spin-degeneracy for more realistic models. 

We describe here the exact DFT solution of the double- 
site Anderson, which corresponds to taking A^ = 1 in Eq. 
(1), and can also be solved analytically. The number basis 
{|n^c.<T >, \nd.a >} of the Fock space is spanned by six- 
teen states, which renders the minimization task of finding 
Q[nc.<j, nd,cT] asymptotically harder. We have found that the 
electron number plane {N^, N^) is split into eight pieces as 
shown in Fig. 4, such that in each piece only a subset of the 
wave-function coefficients is different from zero. As a conse- 
quence, the minimization task has to be performed separately 
for each of those pieces. The Q-functional has again a polyg- 
onal shape. After lengthy algebra, the following expressions 
for the Q-functional in the symmetric case e^ + U/2 = Ec can 
be written: 



a/(1- nc)(l -Ud) + 2x {y^nc,-i - x^ + \/nd.\ - x'^ ) 



+ C/x2, x^ \^|N^ (1 



u 



VTP+WP 



^3 2 



^/(l + me - Ud) (1 + TOd - ?ic) + V(l - "^c - ?^d) (1 - nid - n^) + ^ {U - ^/TPTM^) 



Fi = - 
Fs ^ - 

X 

Fe = - 

Fj ^ - 



v/(l - ric) {l-nd) + 2x ( ^n^,; ~ x^ + ^nd,i -x^)]+Ux^, .t = i W A^^ (l - ^^5^=5) 



v/(nc - 1) {ud - 1) + 2 a; ( v^l - nc,i - x'^ + ^1 - nd,i - x'^ ) + C/ (n^ - 1 + x^), 



VH^J^i-™r-l)(^^r^^nd^^ + v/("-d + TOc - 1) {ric + TOd - 1) 



3-JV 



^C/ - Vt/2 + 64i2) + [/ (nd - 1) 



v/(nc - 1) (7id - 1) + 2 a; (^1 - nc,t - x'^ + ^1 - Ud,^ - ^2 )J +U{nd-l + x^) 



vm+WT- 



Fs = 



^(1 - ne,t) (1 - "d,^) + ^(1 - ne,;) (1 - nd,i) + U (rid - I) 



r 



i21} 



where we have defined a different 



Fa^ 



2\t\ 



for each of the eight a-zones depicted in Fig. 4. We also use 
the site-occupations and moments as rii, rrii — Ui,^ ± n^.^ 
with i = c,d. The full expressions for Fa are shown 
in Appendix B. Simplified expression, valid along the line 
N^ + Ni = 2 are also provided in the appendix. Finally, 
the ground-state energy E^ for given electron numbers N^ is 
found by minimizing Q with the constraints N^ = ric.a + 
Ud^a- To simplify the notation, energies will be measured in 
units of \t\, and the energy origin will be chosen at ec from 
now on. 

We find that Q is spin-degenerate only in regions 1 and 8 of 
Fig. 4, where N is smaller than 1, or bigger than 3. However, 
we find that the spin-degeneracy is lifted if 1 < A^ < 3, be- 
cause here the interplay between kinetic energy and Coulomb 
interactions is more convoluted. The minima of Q and £'" 



occur now along the paramagnetic line M = N- 



T 



N 



i 







regardless of the value of the on-site energy ed, and of U . This 
is shown in Fig. 5(a), where the ground-state energy is plot- 
ted in the [N^, A'^)-plane for the symmetric case and U = A. 
Here the characteristic polygonal shape as well as the pres- 
ence/absence of spin degeneracies in the different regions are 
apparent. The position of the absolute minimum of E'^ along 
the paramagnetic line in contrast does depend on ed and on 
U . For the symmetric case, €d + U /2 = 0, the minimum is 
placed at A^ = 2. Fig. 5(b) shows _E" as a function of N 



along the paramagnetic line for the symmetric case and for 
several values of U which cover the weak-, intermediate- and 
strong-coupling regimes of the model. The chemical potential 
and the energy value of the HOMO are given by the slope of 
these curves. They exhibit the expected discontinuous behav- 
ior at integer values of iViiiMi. Finally, it can be checked that 
this Q functional renders the correct atomic limit by taking 
explicitly t ^- in Eq. (17). As a consequence the exact Q 
functional is free from the delocalization error of mean-field 
theory. The analytic expressions for Q enable to find the exact 
Exchange-correlation potentials for the two-site model {V^^ , 
^cfa)- Notice that these potentials keep the full non-local 
dependence on occupations, because the potential at a given 
site (i, a) depends on all the densities rij^a'- In contrast, it is 
very difficult to determine accurately the non-local terms by 
a numerical solution of this model, or by extending the Bethe 
ansatz LDA approach. We define the exact KS Hamiltonian 
for this double-site Anderson model as 



H 



i—[c,d),(7 cr 

(22) 
This Hamiltonian only has two KS eigenvalues per spin for 
all values of the physical parameters, which are discontinuous 
at integer A^^-values. In other words, the numerical values of 
the KS eigenvalues are constant within each of the eight re- 
gions in Fig. 4, but differ from region to region. We compare 
now the KS eigenvalues with the exact many-body QP spec- 
trum extracted from the poles of the many -body Green's func- 
tion at the impurity's position. Notice again that only integer 




FIG. 5: (Color online) (a) three-dimensional plot of the ground state 
energy as a function of (A'^, Ni) for the symmetric case with f/ = 4. 
(b) Ground state energy along the paramagnetic line in the symmetric 
case for several {/-values. Energies are given in units of \t\. 



electron numbers iV = 0, 1, 2, 3, 4 have a physical meaning. 
For A^ = 1, the system contains a single electron which must 
have either spin up or down. If a ground state wave-function 
with spin up |^i,^ > is chosen, then the spin-up and -down 
Green's functions are different. 



Gd 



Gd 



< *2,t,tl C, 



iti*i.t>p 



Ei^^ — £'2,-t-,'h -\- i 



uj + E^^- Eo + iS 



^i,t 



E 



oj + El.f - E^ +i6 



(23) 



where the summation runs over all spin-0 states with N ~ 2, 
and vE'2,t,t indicate the spin-1 N = 2 state. Similar words can 
be said for A^ ~ 3. Romaniello and coworkers^S have com- 
pared the spectrum of many -body QPs of this model with the 
poles of Green functions evaluated either in the GW approxi- 
mation, or including vertex corrections. They have shown that 
the mean-field static correlation error is amended. However, 
even inclusion of vertex coiTections does not allow to recover 
the QP spectrum in the atomic limit, showing how hard is to 
fully get rid of the delocalization error. 

The exact KS Green's function can be computed using the 
equations-of-motion method giving rise to the following ex- 
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FIG. 6: (Color online) Quasi-particle spectrum of the symmetric 
double-site Anderson model as a function of the electron number A'', 
computed at the thick dots shown in Fig. 4. Green's functions poles 
for f/ = 1 ((al) and (a2)) and [/ = 10 ((bl) and (b2)). The upper 
and lower panels show Gd,t and Gd.j,, respectively. The energies of 
the many-body quasi-particles are shown as black dots whose width 
is proportional to the quasi-particle weight, the Kohn-Sham eigen- 
states are displayed as red dashes. The position of the one-particle 
HOMO level is marked by a thicker dash. Energies are given in units 
of|f|. 



pression 



v: 



xc 



(u 



K5^)(^ 



Cd 






\t\ 



(24) 



This formula must be guided by the results obtained from 
the Lehmann representation at integer N . We compare now 
the poles of the many-body and exact KS Green's functions 
by evaluating V^'-^ at the points in the (iV^, N\^) path shown 
in Fig. 4. This correspond to a paramagnetic solution for 
N = 2 and a spin-up state for N = 1,3. Fig. 6 shows 
the poles of G and G^^ as a function of the electron num- 
ber N for a symmetric case, and for values of U in the weak- 
and strong-coupling regimes. The figure also shows which 
of the KS eigenvalues corresponds to the HOMO level. No- 
tice that the exact many-body and KS spectra closely match 
for values of U not only in the weakly-correlated, but also in 
the strongly-correlated regimes. However, extra many-body 
peaks appear at A^ — 1,2,3, which are not provided by the 
exact KS Hamiltonian. In contrast, the number of many -body 
and KS QPs is the same for A^ = 0,4 because an electron 
added to an empty system or a hole added to a fully occupied 
system can not Coulomb-interact with anything. Occupation 
iV = 2 corresponds to the strongly-correlated Kondo regime 
if U is large, which is the case shown in Fig. 6(b). Here the 
many-body QP spectrum has four poles. These can be clas- 
sified into two sets of peaks placed symmetrically about the 
zero-energy line. The first set is located around ±C//2. The 
two peaks are separated by an energy of order U and corre- 
spond to the upper and lower Hubbard bands. The second 
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set develops into the Kondo resonance for more realistic mod- 
els where Ai is made large. The KS spectrum has only two 
QP, which agree with the two Kondo-like many-body QPs. In 
other words, the KS spectrum shows no trace now of the lower 
and upper Hubbard band precursors. For A^ = 1, the up-spin 
KS spectrum matches the many-body spectrum because the 
many-body phase-space for adding a spin-up electron or hole 
is very limited. However, the many-body phase-space for the 
addition of a spin-down electron is larger, which renders ad- 
ditional many-body spin-down quasi-electron peaks. Similar 
words can be said for A^ = 3, where additional spin-up quasi- 
hole peaks are apparent. Notice in any case that the many- 
body and KS HOMO levels always agree with each other 



V. CONCLUSIONS 

We have presented analytic expressions for the exact den- 
sity functionals of several simple models of strongly cor- 
related electrons, from which we have obtained the exact 
ground-state energy. Those analytic expressions have allowed 
us to write down the full non-local dependence of V'^^ and 
KS Hamiltonians on the occupations. We have computed the 
exact KS eigenvalues and compared them with the true many- 
body QP, as obtained from the poles of the Green's functions. 
We have shown with explicit examples that exact DFT pre- 
serves the quantum nature of electron-electron interactions, as 
opposed to mean-field theory and improvements over it as the 
GW approximation. It is also superior to more sophisticated 
perturbative approximations including vertex corrections. The 



exact functionals do not show any trace of self-interaction, 
static correlation, or delocalization errors. 

We have found that the KS eigenvalues spectrum agrees to 
a large extent, but not fully, with the exact many-body spec- 
trum. This is to say that all KS eigenvalues agree with some 
of the many-body QPs. However, the many-body QP spec- 
trum is richer because the phase-space for addition of quasi- 
electrons or quasi-holes is larger The exact functional only 
warrants the correct position of the HOMO level, while in 
general other KS eigenvalues may or may not agree with the 
exact many-body QP. Remarkably, we have found that the 
KS spectrum most possibly describes the Kondo peak in the 
Kondo regime. However, it is quite plausible that it won't con- 
tain either the lower and upper Hubbard bands, or both. Exact 
DFT has similarities with the Renormalized Perturbation The- 
ory proposed some time ago by Hewsoni^,. The perturbative 
expansion shown in Eq. (2) would possibly describe the full 
many-body spectrum with simple approximations for the self- 
energy. 
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Appendix A: Double-site Hubbard model 



We use the following notation for the hamiltonian of the double-site Hubbard model 



H 






eo 



(ni,o- + n2,a) - ^0 5Z ( ^l.o- ^2,0- + 4,0- C1,<t ) + U ^ fl^. 



tn^,i 



i=l,2 



The expressions for exact density functional are quite similar to those of the double-site Anderson model. 



Fi 

F2 

Fs 



[y^ni,tn2,^ + y/ni^in2,i_ 



v/(l - m +x^- y2) (1 - na - x^ + y^) + (x + y) ( V"i,r - x^ + V'i2,t - V^) + U {x^ + y^) 



^(ni,t - a;2 - z'^) (1 - ni - ^2,^ + :r2 + z'^) + V(n2,; - y^ - z^) (1 - "i,t - "2 + y^ + z'^) 
+ {x + y) (z + ^Jm + 712 - 1 - a;2 - z/ _ ./^M j^ jj (^^i j^ yi^ 



Fi 
F5 



v/(l - ni + a;2 - t/) (l - ^2 - ^2 + y^) + {x + y) (V'^i.i - x^- + ^n^.^ - ?/) +C/(x2+y2) 



^(1 - ni + a;2 - t/) (1 - ^2 - ^2 + y"^) + (2; + y) (^1 - n^.x - iF + \/^ - ?i2,i - x^) 
+[/ (77,1 + 712 - 2 + a;2 + y2) 
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i^fi 



Fv = 



\/("-i,i — 1 + y^ + z^) (2 — ni — 712,^ - y^ - z^) + i/(n2,t - 1 + a;^ + z^) (2 - ni,; — n2 — x^ — z^) 
+ (x + y) fz + v/3 - ni - 712 - x2 - 2/2 _ zA + U (?ii + 712 - 2 + a;2 + y^) 



y(T^7z7+a?2^^^2y|Y^7jJT:^2qr^ + (a; + y) (^l - 77i,^ - y'^ + ^l -n2,^ -a;2) 

+ [/ (?7l + 712 - 2 + x2 + y^) 



F^ 



V/(1 - 77l,t) (1 - «2,t) + V(l - ^hl) (1 - ^i2,i) + C/ (7^1 + 772 - 2) 



where the functionals F„ are defined as 



F„, 



< H > -(?7i +772)eo 



and where Q is found by minimizing Fa with respect to x and y. Along the line iV^ + N^ = 2, the formulae for the F-functional 
can be simplified as follows: 

F = - ( ^rii^i - a;2 + y'l - 772,t - x^j (x + ^/l - m + xA + U {1 - ni + 2x^) 

= - ( y^ni,^ - a;2 + y'l - 772,^ - x^ j (x+y/l-ni+ xA + U {1 - ni + 2x^) 

where the first equation is obeyed if iV^ > 1, iV^ < 1, and vice versa. Q is now obtained by minimizing the above equation with 
respect to x. 



Appendix B: Double-site Anderson model 



The full expressions for Fa are as follows 

F\ = - [ ^ric.t "rf,t + V^c,; n-d,; ] 

F2 = 



^(1 - 77c + a;2 - 2/2) (1 - rtd - x2 + y"^) + {x + y) {y'nZx^xP + \Jnd,^ - t/) 



Uy" 



Fs 
Fi 

F5 

Fe 

Fj 
Fs 



^(7lc,t - x2 - z2) (1 _ ,^^ _ n^.^ + x2 + z2) + v/Ku-y^-^^(l--^^Zt^^^d + y^^^ 

+ (x + y) (z+ yjN - 1 - a;2 



2 9 9 

•" — y — z^ 



C/2/^ 



^(1 - Tic + a;2 - 2/^) (1 - 77rf - 2:2 + y2) + (j: + y) (^77c,^ - x'^ + y/uj^^^) 



Uy' 



v/(l - 77c + x2 - 2/2) (1 - Ttd - x2 + 7/2) + (x + 7/) ( ^1 - "c,i - y^ + ^/l - ^^d,; " a;2) + C/ (n^ - 1 + x^) 



y(77^J^T+y^+^2)(2^77^n^iJ^y^^^^+ v/("-d,t - 1 + a;2 + z2) (2 - ?7c,^ - 77<i - a;2 - z2) 



+ (x + y) (^z + v/3-iV-2:2-y2„^2JJ + [/ (^^ _ i + 3,2) 
v/(l - 77c + a:2 _ y2) (1 - 77rf - a:2 + y2) + {x + y) {^1 - n,,,^ - y^ + ^l-nd,t-x2) ] + U {na - 1 + x^) 



v/(l - 77c,t) (1 - nd,t) + v/(l - nc,i) (1 - nd,i) + U {ud - 1) 



(Bl) 
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Q is again found by minimizing Fa with respect to x and y. Along the line N^ + N^ = 2, the formulae for the i^-functional can 
be simplified, and read as follows: 



F 



rici - x2 + ^/l - Ud,^ - x^j (x+ \J\ - ric + x'^j + U {rid - I + x^) 

j {x + Vl - nc + xA + C/ (nrf - 1 + x^) 



nc,t - a;2 + 4/1 - n^,^ - x"^ 



where the first equation is obeyed if iV^ > 1, iV^ < 1, and vice versa. Q is now obtained by minimizing the above equation with 
respect to x. 



Appendix C: M-slte spinless fermion model 

We show in this appendix the exact DFT solution of the 
double-site spinless fermion model, which corresponds to tak- 
ing Al = 1 and discarding the spin index in Eq. (1), 

H = ecfic + edfid — t{c' d + cf c) + Uhd fic (CI) 

We use a variational wave function of the form 

I > = ao I 0, > + Oc I 1, > + fld I 0, 1 > + fled I 1, 1 > 

(C2) 
to find explicit formulae for the expectation values of H and 
fic, fid as a function of the parameters a^. 



< H > — €cnc + edUd — It cos tp 



\ac\ \ad\ 
D 



<nc> 

<fid> 
D 



I fled I 



D 



\ad\ 



\acd\ 



D 

|aoP + |acP +\ad? + \acd? 



(C3) 



We solve for \ac\, \ad\ in the above equations for nc,d and 
substitute the result back in the equation for < H >. This 



r 



yields the following expressions for < 

-Cd <fld> 



H > -ec < 71c > 



— ZtcosLO i — PI — r — ra 



r, , ^(l-ne)|a^d|2-nd|ao|2 -y/(l-nd)|a^dP-ne|ao |?^ , . 
-'^^('OSip KA'AaoY' ^^> 

where the first and second line apply ifO<Af<lorl< 
N < 2, respectively. The energy functional^ [ric, Ud] is found 
by minimizing the above expression with respect to ao,acd 
and (p. The minimum of the functional happens when Ucd = 
for < A^ < 1, while for 1 < A^ < 2, it is oq which vanishes. 
The resulting functional Q[nc,nd] — ecUc — edfid has the 
following piece-wise shape: 



- 2 \t\ y/ricUd 

- 2 \t\ V(l-nc)(l-nd) + U{n, + nd- 1) (C5) 

where again the first and second line apply if < A^ < 1 or 
1 < A^ < 2, respectively. Q can be easily split into kinetic and 
interacting parts, where both must be defined piece-wise. The 
kinetic term explicitly shows electron-hole symmetry. The in- 
teracting term is non-zero only if A^ > 1, from which a rather 
simple expression for the exact V"^^ can be extracted. 
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